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The conductance of atomic size contacts has a small, random, voltage dependent component 
analogous to conductance fluctuations observed in diffusive wires (UCF). A new effect is observed 
in gold contacts, consisting of a marked suppression of these fluctuations when the conductance of 
the contact is close to integer multiples of the conductance quantum. Using a model based on the 
Landauer-Biittiker formalism we interpret this effect as evidence that the conductance tends to be 
built up from fully transmitted (i.e., saturated) channels plus a single, which is partially transmitted. 

PACS-numbers 73.23.Ad, 72.fO.Fk, 73.40. Jn, 72.15.Lh 

Metallic contacts consisting of only a few atoms can be obtained using scanning tunneling microscopy (STM) or 
mechanically controllable break-junction Jlj techniques. The electrical conductance through such contacts is described 
in terms of electronic wave modes by the Landauer-Biittiker formalism Q . Each of the N modes forms a channel for the 
conductance, with a transmission probability T n between and 1. The total conductance is given by the sum over these 
channels G — X)n=i^n^o 5 with Go = 2e 2 /h the quantum of conductance. By recording histograms of conductance 
values H for contacts of simple metals (Na, Au), a statistical preference was observed for conductances near integer 
values. This statistical preference was interpreted as an indication that transmitted modes in the most probable 
contacts are completely opened (T„ = 1, i.e. saturation of channel transmission), in analogy with the conductance 
quantization observed in 2D electron gas devices [Q. Here, we test this interpretation by performing a new type of 
measurement giving access to the second moment of the distribution of the T„'s. We measure simultaneously, for 
a large number of gold atomic contacts obtained with the MCB technique, the conductance and its derivative with 
respect to the voltage. 
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FIG. 1. Differential conductance dl/dV as a function of bias voltage, measured with a modulation amplitude <0.35 mV, for 
three different contacts with G ~ f .65Go (a), ~ f .02Go (b), and ~ O.88G0 (c). For all three curves the y-scale spans 0.f2 Go- 

The atomic contacts are formed by breaking a gold wire at low temperatures, and then finely adjusting the size 
of the contact between the fresh fracture surfaces using a piezo-electric element [[p. Fig. 1 shows the differential 
conductance, dl jdV measured as a function of bias voltage for three atomic size contacts with different conductance 
values, using a modulation voltage eV <C ksd (with the temperature). For each contact, in order to illustrate 
reproducibility, both the curves for increasing and decreasing bias voltage are given. Measurements such as those of 
Fig. 1 suggest that the fluctuation pattern changes randomly between contact configurations, and that the amplitude 
of the fluctuations is suppressed for conductance values near Go- In order to establish such a relation it is necessary 
to statistically average over a large number of contacts. We do this by measuring the voltage dependence of the 
conductance (dG/dV = d 2 I/dV 2 ) and the conductance itself (G = 31/ dV) by applying a voltage modulation, and 
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measuring the first and second harmonic of the voltage over a resistor in series with the contact. These are recorded 
continuously, while the contact is broken by increasing the voltage, Vp, over the piezo-electric element, producing 
curves as shown in Fig. 2. We use a relatively large modulation amplitude of 20 mV over the contact (at a frequency of 
46 kHz) in order to have sufficient sensitivity and speed of measurement, thus allowing averaging over many different 
contacts. The integration time of the lock-in amplifiers was 10 ms and a reading was taken every 100 ms. Between 
curves the contact was pushed together to a contact with conductance > 20 Go, to ensure that a new contact geometry 
was measured each time. Frequently the contact required manual readjustment and, at these occasions, was pushed 
completely together. All measurements were performed on gold samples of 6N purity, in vacuum at 4.2 K. 
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FIG. 2. Typical example of the simultaneous measurement of voltage dependence of the conductance dG/dV and the 
conductance G, as a function of piezo voltage Vp. The graph includes vertical dotted lines that show that the steps in both 
quantities coincide. Two plateaus have been enlarged and offset to show the tiny steps in the conductance. The elongation of 
the contact is linear with Vp and 10 V corresponds to about 1 nm. 



The conductance in Fig. 2 shows the typical behavior when breaking gold contacts g, which consists of plateaus 
with steps of the order Go and a last plateau close to Go before entering the tunneling regime. The steps and plateaus 
in G correspond with atomic rearrangements and elastic deformation respectively, as the contact is pulled apart and 
finally breaks ||. At each step in G we find corresponding steps in dG/dV . Even tiny steps in G, such as between 7 
and 8 Go, can produce dramatic jumps in dG/dV. 

dG/dV has a random sign and magnitude with a bell-shaped distribution. Fig. 3a shows the standard deviation 
acv = y/((dG/dV) 2 ) — (dG/dV) 2 , as a function of G together with a histogram of conductance values (Fig. 3b) 
determined from 3500 individual curves similar to the one shown in Fig. 2. The squares in Fig. 3a represent the 
standard deviation for a fixed number of 2500 data points each, calculated after sorting the combined data of the 3500 
curves as a function of conductance. Below 0.8 Go where very few points could be measured we use 300 data points 
(circles) . Using fixed numbers of data points in calculating ugv avoids influence of the histogram statistics (Fig. 3b) 
on gov- One clearly observes a very sharp minimum in Fig. 3a at a conductance Go and less pronounced minima near 
2, 3 and even 4 Go- This new observation forms the central result of this paper. Fig 3a shows the combined results 
for three gold samples. The global features reproduce in all three cases, but some sample dependence is observed in 
the shape and height of the maxima. The histogram of conductance values (Fig. 3b) is in accordance with previous 
measurements for gold at low temperature, e.g. M. 

The effect we observe has the same origin as that noted by Maslov et al. J^] in numerical simulations on constrictions 
with defects. The principle can be understood by considering a contact with a single conducting mode having a finite 
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transmission probability T, described by transmission and reflection coefficients t, t' , r and r' (coming in from left 
and right, respectively), with \t'\ 2 = \t\ 2 = T and |r'| 2 = |r| 2 = 1 — T. As illustrated in Fig. 3 (inset), electron waves 
transmitted by the contact with amplitude t, and backscattered to the contact by diffusing paths with amplitude 
a, have a probability amplitude r to be reflected at the contact. This wave interferes with the directly transmitted 
partial wave and modifies the total conductance. A similar contribution comes from the trajectories on the other side 
of the contact. These interference terms will be sensitive to changes in the phase accumulated along the trajectories, 
which is determined by the electron energy and the path length. We can change the energy by the applied voltage, 
giving rise to the fluctuations shown in Fig. 1. Changes in path length of the order of the Fermi wavelength, which is 
the atomic scale, occur at the steps in the conductance in Fig. 2, explaining the correlation with the steps in dG/dV. 
Each time when the contact is opened, and closed again to sufficiently large conductance values, random atomic 
reconfigurations take place, leading to a completely new set of scattering centers. Thus the result presented in Fig. 3 
can be interpreted as the ensemble average over defect configurations. 

In the following paragraphs, we derive an analytical expression for gqy to lowest order in a. In our model, 
the system is divided into a ballistic central constriction connected to diffusive conductors on each side (Fig. 3, 
inset). The central part is described by a transfer matrix t, with elements t nm giving the amplitude for the mode 
n on the left to be transmitted into mode m on the right of the constriction. After diagonalization only a few 
non-zero diagonal elements remain, corresponding to the number of conducting modes at the narrowest part of 
the conductor 0. The t nm are energy dependent, but only on a very large energy scale so that we can ignore 
this in first approximation. For the total transmission of the combined system, in terms of the return amplitudes 
on the left and right hand side of the contact (ai(E) and a r (E) respectively), we obtain the expression t t (E) = 
t;(t /_1 — a r (E)tf~ 1 a[(E) — a r {E)v't'~ 1 — t'~ 1 rai(E))~ 1 t r , with r, r' and t, t' (now in the general multimode case) 
the matrices of reflection and transmission coefficients of the constriction, t/ and t r are the transmission matrices 
through the left and right diffusive regions, respectively. For ksQ <C eV the non-linear conductance can be expressed 
asG = dI/dV, 
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The fluctuations in the conductance are described by 8G = G — (G) (where (G) is the conductance averaged over 
impurity configurations) of which we will consider the voltage dependence dSG/dV. When we take into account 
that scattering processes in the left and right banks are uncorrelated, product terms of a; and a r disappear when we 
average. For the purpose of calculating the small fluctuating part of the conductance we can assume t;tt = t r t£ ~ 1, 
although their deviation from unity will affect (G), which we will address briefly below. Considering first a contact 
with only a single transmitted mode, we obtain an expression for the voltage dependence of the conductance squared, 
averaged over impurity configurations: 
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Products of the form (a(Ei)a* (E 2 )) can be expressed as J Q °° P c i(T)e _1 ' £l_£3 ' T ' ri dT in terms of the classical probability, 
Pci(t) = u f/((1 — cos( , -f))2^/37rkp(DT) 3 / 2 ) to return to the contact after a diffusion time r. We assume the diffusion 
is into a cone of opening angle 7 (Fig. 3, inset), D = vpl e /3 is the diffusion constant and l e — vpT e with r e the elastic 
scattering time [Q. The differentiation of a(eV) in Eq. ([!]) only affects the phase factors (to very good approximation), 
and produces a factor (er/h) 2 under the integral over the diffusion time, r. Further, taking into account that the 
finite modulation amplitude V is the limiting energy scale (ksO, Tl/t^ <C eV with the inelastic scattering time), we 
obtain 
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The T 2 (l — T) dependence results in minima in the amplitude of the voltage dependent fluctuations in the conductance 
at T = or T = 1 and a maximum at T = 2/3. This result can be extended to multiple conducting modes, when we 
assume that the probability to be scattered back to the contact is independent of the mode index, i.e. that defects 
scatter a wave equally into all available modes. The term T 2 (l — T) in Eq. (0) is replaced for the -/V-mode problem 

byEn =1 r„ 2 (i-T n ). 
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FIG. 3. (a) Standard deviation of the voltage dependence of the conductance versus conductance for 3500 curves. The circles 
are the averages for 300 points, and the squares for 2500 points. The solid and dashed curves depict the calculated behavior 
for a single partially open channel and a random distribution over two channels, respectively. The vertical gray lines are the 
corrected integer conductance values (see text), (b) Conductance histogram obtained from the same data set. The peak in 
the conductance histogram at Go extends to 53000 on the y-scale. (Inset) Schematic diagram of the configuration used in the 
analysis. The dark lines with arrows show the paths, which contribute to the conductance fluctuations in lowest order. 

When comparing the experimental data for <tcv with our theoretical model, we need to be aware that the ex- 
perimental data have been sorted according to their conductance value. A given value for G — Go^2T n can be 
constructed in many ways from a choice of transmission values {T n }. The experimental values for <jqv are, therefore, 
an average over impurity configurations and transmission values. Assuming these averages are independent, we can 
compare the data with various choices for the distribution of the transmissions. The dashed curve in Fig. 3a shows 
the behavior of oqv for a random distribution of two T„'s in the interval {0, 1} under the constraint T\ + T 2 = G/Go, 
where the amplitude has been adjusted to fit the data. Alternatively, the full curve shows the behavior for a single 
partially open channel, i.e. in the interval G/Go = {0, 1} there is a single channel, in {1,2} there are two channels 
with one fully open, etc. The latter description works surprisingly well, in particular for the minimum near 1 Go, 
and for the fact that the maxima are all nearly equal. 

Note that the minima in Fig. 3a are found slightly below the integer values. A reduction of the conductance 
with respect to the bare conductance of the contact, Go^T n , results from total probability for back scattering 
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on the same defects which give rise to the fluctuations. We can estimate the correction as the sum over incoming 
channels, n, and their probability to return via any channel, m, J^T n — G/Gq = 2^ n m T n T m {\a nm (E = 0)| 2 ). 

The total return probability (|a nm | 2 ) we approximate by the substitution J°° P c i(t)&t. Thus we expect a correction 
term G = Go(^T n — 2(|a m „| 2 )(^ T n ) 2 ). In Fig. 3 the vertical gray lines indicate the shift below integer values for 
(\a>nm\ 2 ) = 0.005, which is equivalent to a classical series resistance of 130 f2. From this value for (|a nm | 2 ) we obtain an 
estimate for l e = 5 nm, which is consistent with the value obtained from the fluctuation amplitudes discussed below. 

In our experiment we measure the second derivative with a modulation amplitude of 20 mV. This limits the path 
lengths to which we are sensitive to those smaller than ^100 nm (Ly = vph/eV). From the amplitude of the full curve 
in Fig. 3a we obtain an estimate of l e = 5 ± 1 nm, assuming reasonable values for the opening angle 7 of 30°-50° [ pd| . 
This value is consistent with our assumption that d -C l e <C Ly, where d is the contact diameter. The estimate for l e is 
sensitive to the functional form of the factors in front of the T 2 (l — T)-term in Eq. (2), which was not tested in detail. 
Measurements of the dependence on modulation amplitude V are under way. However, the thermopower of atomic 
size gold contacts was recently measured fL2[ and has been found to be determined by the same mechanism, but it 
was measured on an energy scale nearly two orders of magnitude lower. It gives the same estimate of l e = 5 ± 1 nm, 
consistent with the present value. 

Conductance fluctuations |l3| have been observed previously in ballistic contacts with diameters an order of mag- 
nitude larger compared to our contacts, and were measured as a function of both magnetic field and bias voltage p4[ . 
In that work, the quantum suppression of the fluctuations which we report here, is not observable due to the fact that 
many nearly-open channels contribute to the conductance for large contacts. The model introduced by Kozub et al. 
fL5fl to describe the results of Ref. (li) contains only terms due to the interference of two diffusing trajectories, which 
are second order in \a\ 2 . 

The minimum observed at Go in Fig. 3a is very sharp, close to the full suppression of fluctuations predicted for the 
case of a single channel. To describe the small deviation from zero, it is sufficient to assume that there is a second 
channel which is weakly transmitted, T% -C 1, and T\ ~ 1 such that T\ + T% = 1. For this case it is easy to show that 
the value of o~ GV at the minimum is proportional to the average value of T^. Wc obtain (T2) = 0.005, implying that, 
on average only 0.5% of the current is carried by the second channel. For the minima near 2, 3 and 4 Go we obtain 
higher values: 6, 10 and 15%, respectively. 

The well-developed structure observed in acv , with a dependence which closely follows the T 2 (l — T n ) behavior 
of Eq. 2, demonstrates a property of the contacts which we refer to as the saturation of transmission channels: There is 
a strong tendency for the channels contributing to the conductance of atomic-size gold contacts to be fully transmitting, 
with the exception of one, which then carries the remaining fractional conductance. Fig. 3 shows that the positions 
of the minima in (Jgv do not all coincide with those of the maxima in the histogram. This is most pronounced for 
the feature below G = 2 Go- We conclude that the statistically preferred values of the histograms do not necessarily 
correspond with perfect transmission of the bare contact. We propose that the appearance of peaks in the histograms, 
such as the one at 1.75 Go, arises from preferred atomic configurations. 

The concept of the saturation of transmission channels is consistent with recent work, which shows that, for 
monovalent metals, the conductance at G = 1 Go of a single atom is carried by a single mod e [|l6|-[l8[ . Conversely, 
based on the analysis of the subgap structure for superconducting aluminum by Scheer et al. ]17jjl9[ ], which showed 
that typically three channels contribute to the conductance at G = 1 Go, we should expect that aluminum does not 
show a pronounced suppression of conductance fluctuations near integer values. Indeed, preliminary measurements of 
(Jgv on this p- metal, exhibit results for G < Go that are close to a random distribution over 3 transmission channels, 
while the monovalent metals Ag and Cu show behavior similar to Au. 

This work is part of the research program of the "Stichting FOM", which is financially supported by NWO. BL 
and JMvR acknowledge the stimulating support of L.J. de Jongh, and we thank E. Scheer and J.Caro for helpful 
discussions. 
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